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We investigate the dynamics of passive particles in a two-dimensional incompressible open flow 
composed of a fixed topographical point vortex and a background current with a periodic component. 
The tracer dynamics is found to be typically chaotic in a mixing region and regular in far upstream 
and downstream regions of the flow. Chaotic advection of tracers is proven to be of a homoclinic 
nature with transversal intersections of stable and unstable manifolds of the saddle point. In spite 
of simplicity of the flow, chaotic trajectories are very complicated alternatively sticking nearby 
boundaries of the vortex core and islands of regular motion and wandering in the mixing region. 
The boundaries act as dynamical traps for advected particles with a broad distribution of trapping 
times. This implies the appearance of fractal-like scattering function: dependence of the trapping 
time on initial positions of the tracers. It is confirmed numerically by computing a trapping map 
and trapping time distribution which is found to be initially Poissonian with a crossover to a power 
law at the PDF tail. The mechanism of generating the fractal is shown to resemble that of the 
Cantor set with the Hausdorff fractal dimension of the scattering function to be equal to d ~ 1.84. 

PACS numbers: 05.45.Df; 05.45.Pq; 47.52-fj; 47.53 



I. INTRODUCTION 

Large-scale structures such as vortices, background 
and tidal currents are the main constituents in geophys- 
ical flows that strongly influence all the transport and 
mixing phenomena. Transport of heat and mass plays a 
crucial role in dynamics of the ocean and atmosphere. To 
study the transport and mixing processes the Lagrangian 
representation is convenient to adopt. The motion of a 
test particle (a tracer) satisfles the differential equation 



dr 
'dt 



v(r,0, 



(1) 



where v — [x, y) is the position of the particle and v rep- 
resents the incompressible Eulerian velocity fleld. It has 
been recognized that particle trajectories even in simple 
flows can be chaotic in the sense that nearby trajecto- 
ries separate exponentially in time. While in a three- 
dimensional flow it is possible even if the flow is time in- 
dependent 1^ , Lagrangian chaos or chaotic advection in a 
two-dimensional flow may occur only in time-dependent 
flows §111. 

In this paper we study the dynamics of passive particles 
(which take on the velocity of the flow very rapidly and 
do not influence the flow) in a simple two-dimensional 
open flow composed of a fixed point vortex and a back- 
ground current with a periodic component. This model 
is inspired by a very interesting natural phenomena, to- 
pographical vortices over mountains, to be found in the 
ocean and atmosphere 0. In the deep sea, they have 
been found in different regions of the World Ocean with 
the help of buoys of neutral buoyancy § . A laboratory 
prototype of a topographical vortex is a cylindric anticy- 
clonic vortex to be found by Taylor Q in homogeneous 
fluid over an underwater obstacle. In the homogeneous 



ocean the topographical vortices are cylindric whereas 
they are of cone-like form in the stratifled ocean. 

We adopt, of course, an oversimplified model with a 
point vortex in an attempt to catch and quantify the 
main features of the tracer dynamics in the flow from the 
dynamical standpoint. The fixed point vortex is embed- 
ded in a background planar flow composed of steady and 
time-dependent components. If the flow satisfies the in- 
compressibility condition, divv = 0, the evolution equa- 
tion dll) can be written in the Hamiltonian form 



dx , . d'^ 
-=z;.(x, y, t) = -^, 

dy , ,9* 



(2) 



where the streamfunction ^'(a;, y, t) plays the role of a 
Hamiltonian. Thus, the configuration space of an ad- 
vected particle is a Hamiltonian phase space of an asso- 
ciated dynamical system. It means that coordinate-space 
trajectories of tracers coincide with phase-space trajec- 
tories, and more importantly, the configuration space 
contains all the typical structures of Hamiltonian phase 
space including islands of regular motion, islands around 
islands, cantori, stochastic layers, etc. (for a recent re- 
view of Hamiltonian chaos see, for example [|lO|). This 
prominent non-uniformity of the phase space implies 
anomalous transport of tracers, including Levy flights, 
their dynamical trapping and fractal properties, sticki- 
ness to islands boundaries, etc. These phenomena have 
been studied theoretically in the field of chaotic advection 
in a l arg e number of systems (sec, for example, papers 
I, I, y, which represent a small part 

of relevant publications) and observed in laboratory ex- 
periments 11, |g. 

In our model flow, advected particles from the inflow 
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region enter the region where the fixed point vortex is 
located and then wash out to the outfiow region. While 
the tracer trajectories in the inflow and outflow regions, 
where the influence of the vortex is negligibly small, are 
simple as they just follow the streamlines of the back- 
ground current, the trajectories in an influence zone of 
the vortex, called the mixing region in the following, are 
typically chaotic So, the problem of tracer dynam- 
ics in the open flow we address in this paper is, in fact, 
a problem of chaotic scattering |2^ . 

The paper is organized as follows. In Section 2, we 
introduce the model streamfunction and the equations 
of motion of advected particles, consider the integrable 
time-independent version of the model flow and prove an- 
alytically transversal intersections of stable and unstable 
manifolds of the saddle point producing homoclinic chaos 
in the time-dependent system. A numerically simulated 
image of the homoclinic structure and a typical Poincare 
section of the particles trajectories in the mixing region 
provide a general picture of chaotic advection. In Section 
3, we report our main results. We observe and discuss 
the phenomenon of stickiness of tracers trajectories to the 
boundaries of the vortex core and islands of regular mo- 
tion which act as dynamical traps for tracers. We show 
that distribution of the dynamical traps over the space 
of initial tracer's positions is fractal-like, and the trap- 
ping time distribution demonstrates initially exponential 
decay followed by a power-law decay at the probability 
distribution function (PDF) tail with a characteristic ex- 
ponent 7 ~ 2. We compute also dependence of tracer's 
trapping time on their initial positions and find it to be 
a typical fractal scattering function with an uncountable 
number of singularities. The mechanism of generating 
the fractal is shown to resemble the famous Cantor-set 
generation with the Hausdorff fractal dimension of the 
scattering function to be equal to d ~ 1.84. Finally, in 
Section 4, we discuss some possible applications of the 
results. 



II. MODEL FLOW 

We consider a point fixed vortex embedded in a planar 
fiow of an ideal incompressible fluid with a stationary 
and periodic components. The respective dimensionless 
streamfunction 



a particle in the stationary and periodic components of 
the current flowing in the y direction from the south to 
the north. Without perturbation, ^ = 0, the phase por- 
trait of the dynamical system consists of a collection 
of finite and infinite trajectories (streamlines) separated 
by a loop passing through a saddle point with the coor- 
dinates (— 1/e; 0). In the polar coordinates x — p cos 93 
and y ~ p siniy9, the unperturbed equations can be solved 
in quadratures 



e dr 



1 



E-\np 
ep 



-1/2 



dp, 



(5) 



where E ^ ep cos if + liip is the conserved energy. De- 
pending on initial positions and the values of the control 
parameters, particles either are trapped inside the sep- 
aratrix loop and move along closed streamlines or move 
along infinite streamlines outside the loop. In the inte- 
grable system, the stable manifold of the saddle point, 
along which particles move towards the saddle, and its 
unstable manifold, along which they move outward the 
saddle, coincide. Taking into account the separatrix 
value of the energy Eg = —I — \ne, one can easily derive 
the travel time for a particle between two points with the 
coordinates {po, p) nearby the separatrix 



dp' 



Po 



Es + S-lnp' 



(6) 



* = lnv/x2+y2+x(e + ^sinr) = *o+^*i, (3) 



where 5 is a small deviation of the energy from its sepa- 
ratrix value. 

Under perturbation, the stationary saddle point be- 
comes a saddle periodic motion which is represented on 
a Poincare section with the dimensionless period 27r by a 
stationary point. In the extended phase space (x, y, r), 
there exists a stable (unstable) manifold of the saddle 
point with trajectories approaching the saddle periodic 
motion when r — > ±00. Under a typical perturbation, 
stable and unstable branches of the unperturbed separa- 
trix loop split and transverse each other on the Poincare 
section infinitely many times producing a complicated 
homoclinic structure. It can be analytically proved with 
the help of the Poincare-Melnikov integral 



/(«)= / {*o, *i}[2:,(T-a), ?;,(T-a)]dT, (7) 



generates the Lagrange equations of motion of a passive 
particle in the flow Km 



y 



y 



+ y^ 

X 

x^ + y^ 



e + ^ sin r. 



(4) 



where dot denotes differentiation with respect to dimen- 
sionless time T, ^ and e are the normalized velocities of 



where {^'o, ^'i} is the Poisson bracket, Xg and ys are the 
separatrix solutions of the unperturbed problem param- 
eterized by a real number a. With the streamfunction 
(0) we get the function 



/(a) 



OG 00 

sina( J cos r) — cos a ( J drx^sinr), 



(8) 



3 



that obviously has an infinite number of simple zeroes if 
Xs{t) ^ 0. Transversal intersections of stable and un- 
stable manifolds of a hyperbolic point has been analyt- 
ically proven with two-dimensional flows under almost 
arbitrary nontrivial perturbation p3]. A simplified im- 




FIG. 1: Transversal intersections of stable, H^'^\ and un- 
stable, H^~\ manifolds of the saddle point for e — 0.5 and 
^ = 0.01. 



age of such intersections is given in FIG. 1 where a frag- 
ment of the Poincare section with a small perturbation 
^ — 0.01 in the time moments 27rn {n = 0, 1, 2, . . .) for 
a bunch of trajectories, starting from a neighborhood of 
the saddle point, is shown. It is a "seed" of Hamiltonian 
chaos that arises under an arbitrary small perturbation 
around the unperturbed separatrix loop. 

A phase portrait of chaotic advection in our simple flow 
is shown in FIG. 2. The Poincare section of passive par- 
ticle trajectories in FIG. 2a reveals a non-uniformity of 
the phase space that is typical for Hamiltonian systems 
with 3/2 degree of freedom. Fluid around the point vor- 
tex placed at {x = 0, y = 0) forms a coherent structure, 
a vortex core, which is filled with regular, almost elliptic 
orbits. Another large coherent structure is seen to the 
west from the vortex core, it is filled with regular orbits 
as well. A magnification of the northern part of this long 
island is shown in FIG. 2b which demonstrates a chain 
of small islands along its border line. Further magnifi- 
cation of one of these small islands shown in the inset 
of FIG. 2b reveals chains of smaller islands and so on. 
A zone with invariant KAM tori embedded in a stochas- 
tic sea, is called the mixing region. While passive particle 
trajectories just follow streamlines outside the mixing re- 
gion (in the far downstream and upstream regions) they 
can be chaotic inside of it. 
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FIG. 2: Poincare section of the mixing region (a), magnifica- 
tion of the northern part of the long island (b), and magnifi- 
cation of the indicated region between the vortex core and the 
long island in the inset. Parameters are e = 0.5 and ^ = 0.1. 



III. DYNAMICAL TRAPS AND FRACTALS 

Chaotic advection of passive particles in a time- 
dependent open flow is a typical problem of chaotic scat- 
tering. The southern background current with a compar- 
atively small alternating component transports particles 
along the streamlines into the mixing region, where the 
particle trajectories can be very complicated and typi- 
cally chaotic, and then washes out most of them to the 
north where they again follow the streamlines. We show 
in this section that in spite of the simplicity of our flow 
there exists a surprisingly rich variety of different types 
of trajectories inside the mixing region. Particles with 
small variations of initial positions are trapped there for 
trapping times with a broad spectrum extended from a 
particle transit time with the average velocity of the cur- 
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FIG. 3: A chaotic wandering of a passive particle initially 
located at (a;o — —4.358533, j/o ~ —6) in the mixing region. 
The time dependence of the x coordinate (a), the particle's 
trajectory (b) and the particle's locations (black circles) on 
the Poincare section (c) are shown for e = 0.5 and ^ = 0.1. 
The trapping time is T ~ 115. 

rent to infinity. The dynamical trapping cannot be char- 
acterized only by the trapping time T, the ways by which 
advected particles move in the mixing region are very dif- 
ferent. 

Passive particles are placed outside the mixing region 
on the line segment with yo = and different values 
of XQi. We compute the time when the particles reach 
the line y = 6 and positions of the particles, Xi{T) and 
Uiir). The control parameters are chosen to be e = 0.5 
and ^ = 0.1. In FIG. 3 we show the motion of a parti- 
cle with the initial position {xq = —4.358533, uq = —6) 
trapped in the mixing region for the comparatively long 
time T ~ 115. In the upper panel of the figure, the 
dependence of the x coordinate of the particle on time 
T (a) and the particle's trajectory (b) are shown. The 
Poincare plot of the particle positions at the moments of 
time T = 2?™ (black circles) is shown in FIG. 3c. The 
unperturbed separatrix and the Poincare section for a 
collection of particles initially placed inside the mixing 
region are shown for reference. The particle moves quite 
uniformly without any preference to visit special zones in 
the mixing region. 

The second particle to be placed initially close to the 
first one at the point (a;o = —4.358034, yo = —6) demon- 
strates cardinally different behavior shown in FIG. 4. 



It spends in the mixing region much more time, T ~ 
1581, than the first particle does. The figure clearly 
demonstrates the phenomenon of stickiness well-known 
in chaotic Hamiltonian dynamics [|o[ In con- 

trary to the first particle, the second one spends most of 
time in the mixing region nearby the boundaries of the 
vortex core and the long island. Motion of the particle 
in these areas is almost regular with zero maximal Lya- 
punov exponent. In Hamiltonian systems, there may ex- 
ist partially broken invariant curves (cantori) in a vicinity 
of island boundaries which prevent particles to quit the 
area and act as a kind of permeable barriers. Both these 
reasons result in stickiness of the islands boundaries that 
is seen in FIG. 4. 
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FIG. 4: Stickiness of a passive particle initially located at 
(xo = —4.358034, j/o = —6) on the boundaries of the vortex 
core and the long island. The time dependence of the x co- 
ordinate (a), the particle's trajectory (b) and the particle's 
locations (black circles) on the Poincare section (c) are shown 
for e — 0.5 and ^ = 0.1. The trapping time T ~ 1581 is very 
long. 

To study the trapping time distribution we have car- 
ried out a number of numerical experiments. A trapping 
map referred to initial positions of = 22801 passive 
tracers, distributed homogeneously on a rectangular grid 
{xq e [—5, 0]; yo G [—6, 1]), is shown in FIG. 5. Color 
modulates values of the trapping time T. The inset shows 
the respective zone of this map magnified by the factor 
of 30 on the x axis and factor of 17 on the y axis. A 
patchiness exists for all the scales revealing a fractal-like 
structure of the underlying phase space. Tracers started 
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FIG. 5: Trapping map showing dependence of the trapping times in the mixing region T on initial positions of A'^ = 22801 
tracers distributed on a rectangular grid in the upstream region. The inset shows a zoom of the fragment of the map. 



initially in the areas of the upstream region, correspond- 
ing to black patches on the map, will be trapped in the 
mixing region for longer times. The inset demonstrates 
that there exist tracers with very close initial positions 
whose trapping times differ largely. To quantify this set 
of dynamical traps with a wide range of trapping times 
we compute the respective distribution function N(T) for 
N = 110011 tracers (positioned initially on a rectangu- 
lar grid {xo e [-4.65, -4.35]; yo £ [-6.05, -5.95]). The 
double logarithmic plot in FIG. 6 demonstrates initially 
exponential decay followed by a power-law decay at the 
PDF tail, N{T) ~ T^"', with the characteristic exponent 
7 ~ 2. While a Poissonian distribution is expected for 
strong chaotic mixing a power-law dependence indicates 
the presence of dynamical traps in the flow. Particles 
trapped for longer times give an enhanced contribution 
to the PDF tail. Different choices of initial positions of 
the tracers visiting the mixing region do not change sig- 
nificantly the trapping-time distribution. 

The well-known manifestation of chaotic scattering 
is fractal-like scattering functions with an uncountable 
number of singularities [|l6|, |l^, ^ . We take the trapping 
time to be the scattering function, T{xq), and compute 
it carefully for N = lO** particles distributed initially on 
the line segment with = —6 and xq S [—4.65, —4.35]. 
The results in FIG. 7 demonstrate a typical scattering 
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FIG. 6: Trapping time distribution with A'^ = 110011 tracers 
shown on log — log axes. 

function r(a;o) with an uncountable number of singu- 
larities which are unresolved in principle. The inset in 
FIG. 7a shows a zoom by the factor of 20 of one of these 
singularity zones. Successive magnifications confirm a 
self-similarity of the function with increasing values of 
the trapping times T. FIG. 4 gives an example of a tra- 
jectory with very long trapping time T ~ 1581. Remind, 
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FIG. 7: Fractal dependence of the trapping time T on initial particles positions xq with the inset showing a 20-fold magnification 
of one of the singularity zones (a). Mechanism of generating the fractal with magnification of a small segment corresponding 
to the inset in the top panel (b) . 



please, that we measure time in units of the period of the 
alternating current. The trapping map in FIG. 5 provides 
a two-dimensional image of the fractal structure of the 
scattering function T{xq, ya). As wc have shown above, 
sticky boundaries of islands of regular motion, embedded 
in a stochastic sea, act as a kind of dynamical traps pro- 
viding all the spectrum of values of the trapping times 
up to infinity. Namely these dynamical traps are the ul- 
timate reason for appearing fractals in our simple model 
flow. To give an insight into a mechanism of generat- 
ing the fractal we plot in FIG. 7b segments of the initial 
string with 10^ tracers which arc trapped in the mixing 
region after n rotations around the fixed point vortex. 
After each rotation, a portion of tracers is washed out in 
the downstream region with y ^ 6. This process resem- 
bles the mechanism of generating the famous Cantor set 
but is more complicated. 

In conclusion, we compute the length of the curve 



T{xo) shown in FIG. 7a 

N 

L{l/N) = Y,\T,+r-Ti\, (9) 

i=l 

as a fimction of the inverse number of tracers, , dis- 
tributed initially on the same line segment with yo = ~6 
and xo G [—4.75, —4.25] in the upstream region of the 
flow. It is a decreasing (in average) function with a neg- 
ative slope whose double lagarithmic plot is presented 
in FIG. 8. The Hausdorff dimension of the fractal-like 
scattering function is computed to be d ~ 1.84. 



IV. CONCLUSION AND DISCUSSION 

We have shown that chaotic advection of passive par- 
ticles by the open flow, composed of a fixed point vortex 



7 



L 




FIG. 8: Length of the fractal curve T{xo) versus the inverse 
number of tracers A'^"^ distributed initially on the same line 
segment in the upstream region. The double logarithmic plot. 

and a background current with a time-dependent com- 
ponent, has a fractal nature. Appearance of fractals and 
dynamical traps in geophysical flows with topographical 
vortices should fluence strongly transport and mixing of 
heat and mass. Instead of homogeneous mixing, one ex- 
pects a highly-structural, hierarchical fractal distribution 
of advected particles that is of great interest in oceanog- 
raphy and atmospheric sciences. It should be emphasized 
that the fractals and dynamical traps are not specific 
for our particularly simple model. They are typical in 
low-dimensional Hamiltonian systems and should exist 
in more realistic models with elliptic topographical vor- 
tices and different boundary conditions considered, for 
example, in p^ . 



If tracers are not passive but active (chemical reagents 
or biological species), their activity is superimposed on 
chaotic mixing in open background flows. Typically, 
there exists a local coupling between a spatial charac- 
ter of the activity and the spatial character of the fluid 
with coherent structures, dynamical traps and fractals. 
Say in chemistry, molecules will react if they are in a 
close contact with each other. So, the imperfect mixing 
should largely increases the chemical or biological activ- 
ity in some regions of the flow and suppress it in other 
regions. In marine biology, the paradoxical observational 
fact of coexisting thousands of phytoplankton species in 
spite of a limited number of resources for which they 
compete has been known for a long time ||2^. Chaotic 
mixing of active particles in open nonstationary flows in 
the ocean with small-scale spatial inhomogeneities and 
fractal patterns can provide a natural explanation of the 
plankton paradox . Homogeneous and complete mix- 
ing leads to an equilibrium state in the system of com- 
peting species resulting in the survival of the most fitted 
species, whereas chaotic mixing generates, typically, a 
persistent nonequilibrium state in the system providing 
the coexistence of much more species than the number of 
ecological niches they occupy. 
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